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Abstract 

Dynamics of a discrete polymer in time-dependent external potentials is studied with the master 
equation approach. We consider both stochastic and deterministic switching mechanisms for the 
potential states and give the essential equations for computing the stationary state properties of 
molecules with internal structure in time-dependent periodic potentials on a lattice. As an exam- 
ple, we consider standard and modified Rubinstein-Duke polymers and calculate their mean drift 
and effective diffusion coefficient in the two-state non-symmetric flashing potential and symmet- 
ric traveling potential. Rich non-linear behavior of these observables is found. By varying the 
polymer length, we find current inversions caused by the rebound effect that is only present for 
molecules with internal structure. These results depend strongly on the polymer type. We also 
notice increased transport coherence for longer polymers. 
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I. INTRODUCTION 



There has been considerable progress in the research of Brownian motors during the last 
decade (see e.g. jl, 2, 3|). Starting with the simple pointlike Brownian particles with time- 



dependent driving forces, research hasjsxjmndec 
interacting Brownian particles (e.g 



has expant 



towards the more comp 



9| ) and polymers 



10 



11 



ex objects such as 



121 ] . In this paper, 



we study polymer motion with discrete lattice models, which allows us to consider different 
kinds of microscopic polymer dynamics in detail. Aside from being a purely theoretical 
branch of study, analysis of simplified discrete non-equilibrium particle models has became 
an important tool for studying biologically inspired Brownian motor systems (e.g. [13]' 
discrete models have been applied widely to single particle ratchet problems (e.g. [14 



16 



171]) . We expand this picture by considering a generalized Rubinstein-Duke model (RD 



model) 18|, ll9| for polymer motion in discrete time- dependent potentials. An interesting 



question is, what kind of dynamics lies beyond simple pointlike particles and how one can 
calculate its properties such as the effective diffusion coefficient and the drift. Although 
there are plenty of studies concerning the behavior of the RD polymer in zero (or uniform) 



field (e.g. 



20 



21 



in this context 



]), only recently a ratchet mechanism (tilting ratchet) has been considered 



22|. 



Especially because of the high complexity of Brownian motors with internal structure, 
most studies of these systems have applied the Monte Carlo method. However, since the 
ratchet systems are quite sensitive to the values of parameters, and drifts generated by pure 
ratchet mechanism are usually very small, Monte Carlo simulations tend to be very time- 
consuming and inaccurate. In this paper we study these systems with the master equation 
approach. The results obtained this way are accurate enough to reveal the details of the 
dynamics. 

The purpose of this paper is to give a hands-on example of how one applies the mas- 
ter equation method to systems involving time-dependent periodic potentials and complex 
molecules by using a modified RD model polymer as a prototype of such molecules. We 
perform calculations for short linear polymers in a non-equilibrium environment generated 
by flashing and traveling ratchets. To test the significance of the polymer type, reptating or 
not, we compare the motion of the RD polymer with the dynamics of a modified version of 
the RD polymer with less constrained microscopic movement. 
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The paper is organized as follows: In Section II we expand the modified RD polymer 
model to periodic time-dependent potentials and give equations for the calculation of the 
drift and diffusion coefficients, in Section III we present results of the calculations for short 
polymers and finally in Section IV we give our conclusions and discuss the implications to 
applications. 



II. THE MODEL AND METHODS 



The RD model was originally developed to model the random motion of a flexible polymer 
in a confined medium with static obstacles (e.g. pores in gel) that the polymer must bypass, 
therefore causing the polymer reptation. By assuming that the network of obstacles can 
be modeled (on average) by a lattice-like structure, that the correlation length between the 
polymer segments is smaller than the distance between the obstacles, and that only the 
polymer heads are able to move into previously unoccupied cells (lattice sites), the problem 
can be discretized to a staple particle hopping mode, Q. Soon after the original model 
was expanded [19( to be suitable for external potentials (e.g . static field), pure theoretical 
research of the model started to flourish such as in Refs. 20|, 123 ]. 

Technically the RD model is a spin-1 chain with special kind of nearest neighbor inter- 
actions between the particles (reptons). By assuming that the reptons experience random 
"pushes" by the environment modelled with a continuous time Markov process with expo- 
nentially distributed waiting times, we can construct the stochastic generator of the system. 

In order to study the effect of the intrinsic transition rules of the polymer in time- 
dependent periodic potentials on long-time dynamics, we will compare the results of the RD 
model to the results of a non-reptating polymer which allow the breaking of the reptation 
tube. In this paper we call this extended model the free-motion model (FM model). In 
Fig. 1 there is an illustration of an example configuration of a six repton polymer with 
arrows indicating all allowed moves for both RD and FM models (see also Ref. 24j). All 



repton transitions are between nearest neighbor lattice sites 



24 



25 



only . Similar extensions have 



26|. 



been studied previously in a different context in Refs. 

As an environment for the polymers, we assume a discrete periodic potential V(x) such 
that 

V(x + L) = V(x). 
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FIG. 1: (Color online) Illustration of the allowed transitions in RD and FM models for a six repton 
polymer in one of its configurations. The moves described by the blue arrows are only allowed 
for the FM polymers and those by the black arrows for both polymer types. The letters a, b, c 
represent the operators corresponding to the moves and are defined later in the text. 

To make contact with Kramers rate theory (see e.g. review [^j) and the previous work 
related to discrete ratchets [16j, we define the transition rate from state i to j by 

7 ^ = aexp (p(V(i)-V(j))), 

and choose a = (3 = 1 to define the time and energy scales along with the lattice constant 1 
to define the spatial length scale. We shall next define the time-evolution operators for the 
RD and FM models (readers not interested in the formal development may skip the rest of 
this section). 

The mathematical model for the polymer, which contains the RD model as a special 
case but also allows breaking of the reptation tube if wanted, is constructed as follows 

n 

(see e.g. |28(). Within the most compact, the inner coordinate representation, every bond 
between reptons can be in three states; up (state A), down (state B) or flat (state 0). In 
Fig. 1, reptons 1, 2 and 3 are in state 0, repton 4 in state A, and reptons 5 and 6 in state B. 
An iV-repton polymer has N — 1 bonds. The state corresponding to polymer configuration 
y is thus given by a 3 Ar_1 -dimensional state vector \^f y ). 

The non-zero elements of the local creation and annihilation operators defining the dy- 
namics of the bonds are 

Mi.i = M 2 , 2 = M 3 ,3 = 1 

Wa,i= Hi, 2 = [ & ] 2 ,3= [ &t k 2 = 1 - 

The operators a and b produce changes in the local bond configuration as indicated in Fig. 1. 
To extend the model to include a periodic potential V, we must add an additional state. 
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One repton is chosen as a marker repton that keeps track of the polymer position within the 
potential. The transition rates of a single repton now depend on the position of the marker 
repton and all other bonds separating it from the marker. Either of the head reptons is the 
most convenient choice for the marker repton, hence we choose here the repton labeled 1 
(see Fig. 1). The dimension of the marker state is L, so the dimension of the total system 
of equations becomes L x S^" 1 . 
By denoting 

L(i)=exp(-V(i + 1) + V(i)) 
R{i)=exp(-V{i-1) + V(i)), 

indicating transitions to left and right (corresponding down and up in Fig. 1), the non-zero 
matrix elements for the marker state and transition operators are 

[cr],_ M = 1 for ^ 1, h + ] I+1>| = lforZ^L 

= H]i,l = Mm = 

where 1 < I < L. The state of the polymer now has the form 

(marker repton) <g> (polymer configuration) = |^) ® \^ y ), 

where is the marker repton state vector with dimension L. The stochastic generator of 
the polymer model in the L-periodic potential thus becomes 

L 

i=i 

where the operator A applies to bond 1 and the marker repton, M applies to bulk reptons 
and B applies to bond N — 1. The explicit forms of these operators are given in Appendix 
A. 

A. Time-dependent potentials 

The time-dependence of the environment can break the detailed balance and may result 
in a directed drift. We assume that the switching between the distinct environments is 
independent of the polymer state in the potential i. e. there is no feedback from the polymer. 



N-2 



i,y,l 



(1) 



i=l 
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The switching mechanism between the potentials can be either stochastic or deterministic. 
The stochastic Markovian switching allows us to evaluate the stationary state directly by 
solving an eigenvalue problem. This is the most widely used way of studying Brownian 
motors and similar systems. With deterministic switching, we must numerically integrate 
to get the periodically stationary state. 

Due to the time- and position-dependent transition rates, extra care must be taken to 
numerically study the process accurately. For example, one should not use the standard 
discrete-time Monte Carlo simulation method that has been widely used in various RD 
model studies. It does not produce correct results for our models. Instead one should handle 
the master equation directly by means of numerical integration or use the continuous-time 
Monte Carlo method. In general, transport of particles in time-dependent potentials is a 
hard problem to solve exactly. Even for a single particle in a periodic potential the general 
solution is not known. The solution for stationary potentials, however, is available (pioneered 
by Derrida |2s|). 

First assume Markovian switching between the potentials. We must include an additional 
state that keeps track of the current potential 

(potential state) ® (marker repton) ® (polymer configuration) = \ty s ) ® <g> \^ y ), 

where |\l/ s ) is the state vector of the potential with dimension S i.e. the number of different 
potentials. Since there is no feedback mechanism, adding this new state is straightforward. 
The non-zero state and transition operator elements for the potential state are 



= 1 for 1 < i < S 



h, 



= l n s s 
J 1,5 



where 1 < s < S. By defining the operator h like this, we consider only cyclic transitions 
between the potentials (i.e. 1 — > 2 S 1 —►...) to preserve the analogy with the 

deterministically switching potentials. The stochastic generator becomes 

s 

H = J2[H s + T- 1 (n s -h s ^ , 

s=l 

where if s 's are formed by extending all elements of the operator in Eq. ([1]) with their 

corresponding potential state s (e.g ay — > n s aij) and Ti,T 2 , . . . ,T$ are the mean life- 
times of the potentials. 
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With deterministic switching, the stochastic generator is given by 

H x , t G [0, Ti) 



H(t) 



H 2 , tE[T 1 ,T 1 +T 2 ) 



t e 



where H s is the operator of the type ([!]) in the potential s and T = Y^i=i Ti is the time- 
period and H(t + T) = H(t). In this case there exists a T-periodic stationary solution. 
Once H is given, the time-evolution of the system is governed by the master equation 
dq(t)/dt = H(t)q(t), where q(t) is the probability vector. Since if s 's are generally non- 
symmetric, q(t) usually has an oscillating behavior. 

B. Drift and diffusion 



We are interested in the drift and diffusion of the center of mass of the polymer. The 
velocity and the diffusion coefficient can be defined as 

d 



D, 



eff 



lim —{x C u(t)) 
t^oo at 

1 lim ^ ((x CM (t) 2 ) - (x CM (t)) 2 ) 

2 t^oo at v 



where xcm is the center of mass of the polymer. Here v and -D e ff could also be defined for 
single reptons instead of the center of mass and this local approach naturally leads to the 
same longtime values. 

From the previous we define the Peclet number 

\v£\ 



Pe 



where we choose the length scale I = 1. Since our polymer is simply composed of sev- 
eral neighbor-hopping random walkers, we can generalize the formalism of Ref. 17| (which 



generalizes the ideas of Ref. 



29) ). First define 



n=—oo 
oo 

S ?/00 = i l + nL )Pn,l,y'{t) - (x C u{t))q y {t), 



(2) 



(3) 
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where p n ,i,y is the probability to find the marker-repton in the position I + nL with the 
polymer inner configuration y' and the re-defined state y includes both the marker-repton 
position (I) and the inner configuration (y f ) within the L-periodic potential. Assume that 
the stochastic generator H of the total system is defined by the rates F^j from state i to j. 
It can be shown by using the definitions above, by taking the time-derivatives and using the 
master equation (see e.g. {30] for a similar calculation) that 

v(t)=h^2 (R V ont - iLt) Qy(t), 
y 

where 

-^out — / , L y,i -^out ~~ / j L y,i 

i i 

with arrows indicating the direction (right or left) of those repton transitions that lead from 
the state y to states i, neglecting all the rest. Since this expression is for the center of mass, 
h = 1/N is chosen as the new lattice constant. 
Similarly we get 

y 

+h^2(R y ont -L y ont ) s y (t). 
y 

The evolution equations for q y (t) and s y (t) can be found by differentiating ([2]) and ([3]) in 
time and using the master equation once more. We arrive at 

dj ^ = - PC + IL) %(t) + R*Mt)) + LfMt)) 

dSy(t) -jD ^ _ w -j^y / _ / ±\\ , TV 



(Jf - (R y out + L y out ) s y (t) + Rl(s(t)) + L y m (s(t)) 
h[L»Mt))-K(Q(t))] -v(t)q y (t), 



where 



i i 
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Note that all transitions are assumed to be between nearest neighbor lattice sites only. 
Otherwise transitions of certain length should be collected in their own sets according to 
their hopping distances (< L), which would appear as coefficients of additional sum-terms 
in the equations. In the matrix form 

dq(t) 

dt 
ds(t) 

dt 

where H siRn has the structure 



H(t)q(t) (4) 
H(t)s(t) - hH sign (t)q(t) - v(t)q(t), (5) 



[H sign ] iyj = [H} iyj for all right transitions 
[i/sign]jj — — for all left transitions 
[i/sign]jj — for all other transitions i,j . 

This operator is easily built while building the stochastic generator itself. Since v(t) in 
Eq. (jSj) is governed by Eq. (jlj), these systems must be solved simultaneously. See also 



Refs. m 



32| where similar approach has been applied to find the drift and the effective 



diffusion coefficient for complex molecules. 
1. Time-independent stationary states 

When H is time-independent, we can take the limit t — > oo and define the steady-state 
parameters as 

Q y = lim q y (t) S y = lim s y (t). 

t — >oo t — *oo 

By using these we get well-defined stationary values 

v = h Yl K>t - ^Wt) Qy 
y 

D cS = C^out + -^out) Qy 

y 

+hJ2(Xut-lZut)Sy (6) 

y 

for the velocity and the effective diffusion coefficient. Now Q y s and S y s are found by solving 
the equations 

HQ = HS = hH sign Q + vQ. (7) 
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So far equations like these have been solved exactly only for a single particle on a periodic 



lattice. The first solution was given in Ref. 



29] for the nearest neighbor hopping particle 



with arbitrary transition rates. This has been later extended e.g. for parallel one-dimensional 
lattices in Ref. 3(3]. However, for more complex systems (like RD polymers), solutions cannot 
be found by exact methods and numerics must be applied. The structure of H also raises 
some issues. Since the determinant of H is always zero, mathematically there is no unique 
solution for the non-homogeneous linear set of equations in ([7j). This can be easily seen by 
using the fact that an ergodic stochastic system always has a non-trivial stationary state, 
therefore by the rank-nullity theorem we have Kank(H) = Dim(iJ) — 1, meaning that we 
have one free parameter and all solutions are separated by a constant (i.e. S y is a solution 
<=>■ (S y + constant) is a solution). However, since we also have conditions 

y y 
which can be derived from the definitions of s y (t) and q y (t), there indeed exist unique 

solutions for S and Q (which is of course also required on physical grounds). 



Eq. ([6]) is a generalization of the result derived in Ref. |20j . This can be seen by considering 
the case L = 1 without external potentials (i.e. y's are simply inner configurations, v = 
and Q := Q y = 3~ N+1 ), so that with a y := —2S y /hQ we have (for lattice constant 1) 



D eS = -Qj2 



( — y — y \ ^ { — y — y \ 

(-^out + Ajut) + J^q (^out ~ Aut) Sy 



2iV 2 3 



v 

»JV-1 



Here we used the fact that, in this case, every state y has a weight 1/3 and r y /lyS can 
be in terp reted as the "number of arrows" for right/left transitions out from the state y. In 
Ref. }20l | . symmetry properties (reflections) of polymer configurations were used to find a 
unique solution for a y 's, but this is not possible when external potentials are present and 
the problem is non-symmetric. However, numerical linear algebra tools can be used to find 
the solution. 

2. Time- dependent stationary states 

When H is time-dependent, we must integrate equations (@J and (jSJ) in time until the 
system arrives at the periodically stationary state (with period T). The mean velocity and 
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the diffusion coefficient are determined by 



v = lim — / v(s) ds 



t—>oo T 



t-T 
t 



D cS = lim ^ / D eS (s) ds. 



t-T 



In practice, these are calculated by integrating in time long enough so that results have 
converged. 



C. Fast and slow switching regimes 

When the switching times of the potential are close to the characteristic timescales of 
the system (i.e. relaxation times), the behavior depends heavily on the switching type and 
lifetimes of the states. However, when the potential changes very rarely or extremely fast, 
the system becomes independent of the switching type and even of the relative life-times of 
the states. 

First assume that the total mean switching period T — > such that Tj > for all mean 



lifetimes of the potentia 
potential ("mean-field" 



s (T < i < S). In this case particles experience an effective average 
361 ]) and the transition rates become 



r!? = (8) 

k 

where r* „• are the transition rates of the stochastic generator of type (CQ) in the potential 
k and Xk = Xfc/T's are weight factors determined by the mean life-times of the potentials. 
This leads to a mean- field stochastic generator with dimension L x S^ -1 . This approach was 
used in Ref. 33J to solve exactly the single particle dynamics in two arbitrary alternating 
periodic potentials. Although this mean-field limit is mathematically well defined, from the 
physical point of view it's artificial since real- world systems have inertia, and changing the 
potential state takes some finite time (e.g. charge re-distribution to build up an electric 
field). So the velocity always goes to zero in the fast switching limit. 

Now assume that Tj ^> Tj for all 1 < % < S where r» is the longest relaxation time of 
the system in the potential i. This means that the system always converges close to the 
stationary state in the current potential before the potential is switched to the next one. By 
the model assumptions, drift is always zero at the stationary state in all potentials. Let dju 
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denote the mean travel distance of the molecule center of mass within the potential j using 
the stationary state of the potential i as an initial state and then letting the system fully 
relax [3?]]. Summing over all dj^s gives the total expected distance within one time-period 
T, and by assuming cyclic switching of the states, we define 

5 

d = X) ^K+i- ( 9 ) 

i=i 

The sign of d determines the drift direction in the large T limit and the asymptotic drift 
thus becomes v = d/T. 

That internal molecular states may have strong influence on the dynamics can be already 
seen in the slow switching regime. Letting the molecule first find its equilibrium in some 
non-flat potential and then turning the potential off may indeed result in directed motion of 
the molecule after the switching, due to internal relaxation, whereas a single particle would 
be immobile in the mean. These rebounds might be dominating and define the sign of d. 



III. RESULTS 



A. Choice of the potentials 

We have numerically analyzed RD and FM models with the polymer length of N = 1...11 
reptons and with two potentials (S = 2) and stochastic switching. All calculations were done 
with MATLAB. We used the standard Runge-Kutta 4 method to integrate fl3J) in time and 
a trapezoid method to calculate the resulting integral in Q. The Arnoldi and BiConjugate 
gradient stabilized methods were used to solve homogeneous and non-homogeneous systems 
in ©. 

We consider two basic potential types: flashing and traveling ratchets. The first type 
is the most general non-symmetric potential that has been extensively used in studies of 
Brownian motors and the latter one is a generic example of asymmetrically placed symmetric 
potentials and has been recently used with single particle models [17|, |34|. We consider the 
simplest case L = 3, which is the smallest possible length that can form both of these 
potentials with the ratchet effect. See Fig. 2 for sketches of these potentials. A positive drift 
sign indicates motion in the increasing lattice-site index direction. Because of symmetries, 
the next choice would be L = 5, but this choice would also need longer polymer lengths 
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FIG. 2: Sketch of the flashing non-symmetric ratchet (left columns) and traveling symmetric ratchet 
(right columns) for L = 3 (two period lengths shown). 

(N ^> 11) than we can efficiently handle. We require that the polymer must be able to 
cover several potential periods when fully extended. We also set Vm ax = 1/2, which we 



found to give interesting results while also being computationally feasible [38|. The results 
concerning the general behavior and drift inversion do not significantly depend on the choice 
of Vmax- With these parameter choices the relaxation times r for iV = 3. ..11 fall between 
In (r) = 0.8, ...,4.8 for RD polymers and In (r) = 0.3, 3.4 for FM polymers in all potentials 
studied here. Below we let N and T vary. By the limit T — > we mean going to the mean- 
field stochastic generator with the rates given by (jHJ). Overviews of the dynamics of polymers 
of lengths N = 1,3, 5, 7, 9 are given in Figs. 3 and 5, while Figs. 4 and 6 provide more detail 
for iV = 1...11. We are especially interested in the current inversions and the general effects 
of the polymer size. 



B. Flashing ratchet potential 



Let us first define the time-period T = T on + T Q s and the symmetry parameter x = T Q g/T, 
where T on / g are the corresponding mean life-times of the potentials (see Fig. 2). In previous 



studies (e.g. [15 



16]), only symmetric flashing x = 1/2 was considered. This results in zero 



drift for T — > 0, which also happens in all real systems (for all x). However, with x ^ 1/2, 
this does not happen for the models considered here. The drift changes its sign as a function 
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FIG. 3: (Color online) Drift and diffusion of the RD and FM polymers in the flashing ratchet 
{L = 3, Vmax = 1/2)- (a-c): drift as a function of the symmetry parameter x = T g/T with the 
total flashing period T — ► (a), T = exp (3) (b) and T = exp (5) (c) with N = 1 (solid black), 
N = 3 (dash), N = 5 (dot), N = 7 (dash-dot), N = 9 (solid blue), (d): effective diffusion 
coefficient of N = 9 with T — > (solid), T = exp (3) (dash) and T = exp (5) (dot). 

of x and the point of this sign change in x depends on T. This is shown in Fig. 3, where we 
have plotted v as a function of x with three different T's (figures (a)-(c)) that represent the 
general behavior in different scales of T. The drift in the positive direction (generated by 
the short slope) arises when the ratchet is switched on for such a short time that the larger 
rate of the short slope wins the smaller rate of the longer slope (see Fig. 2). Therefore, for 
increasing T, the ratio x must get smaller to retain the dominance of the shorter potential 
slope, and finally x goes to zero at T — > oo. 

When we add more reptons, the overall shape of the v curves remains very similar with 
small T's. However, a clear effect of the polymer length and internal mechanisms can be 
seen with the long-time period T = exp (5), where the drift curve of the N = 9 RD polymer 
turns positive for x ~ 1/2. The velocities of the FM polymers remain on the negative side 
and no change in their general drift behavior can be seen as the parameters N and T are 
varied. In Fig. 3 (d) we have plotted an example of _D e g behavior of the N = 9 polymers with 
corresponding T's of the figures (a)-(c). Here x = 1 simply gives the diffusion coefficient 
of the free polymers, and the diffusion constant in the static potential (at x = 0) is always 
smaller. As can be seen, the effect of T and x on the diffusion is quite small in general. 
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Next we fix values x = 1/4, 1/2, 3/4 and examine the T dependence of the drift and 
the Peclet number in detail. The results in Fig. 4 reveal a complex behavior of the drift. 
The overall form of the curves is as expected: the drift and the Peclet number have some 
(local) maxima around In (T) « 0. For small T, the single particle remains the fastest in all 
cases excluding x = 1/2 for FM polymers, where it is the slowest one. However, as T gets 
larger, longer polymers eventually become faster, which is caused by their longer relaxation 
time (short polymers have already reached their stationary state). This can be clearly seen 
from Figs. 4 (b) and (c), but it also takes place in figure (a) to some extent. Similar behavior 



of coupled partic 
reported in Ref. 



lies being faster than single ones and also having drift inversions were also 
5j. Although the relaxation times are quite different (see Sec. Ill A), the 
maxima of the drift fall close to In (T) « 1 for all polymer lengths and the position of the 
maximum Peclet number is almost constant. The drift sign change, already seen in Fig. 3, 
is present in Fig. 4 (c). 

The behavior of the Peclet number is very clear and similar in every case in Fig. 4: the 
larger the polymer, the larger the Peclet number. Thus the transport of longer polymers 
is more coherent than of shorter ones. Similar behavior was found in a continuum model 
consisting of elastically coupled Brownian particles [4|. By comparing the values of the 
Peclet number between polymer types, we see no significant differences between the curves. 
There is a slight difference for large values of T, where the Peclet number remains larger for 
FM polymers. This holds with every choice of parameters, excluding the possible current 
inversion points (e.g. the interval In (T) = — 1...0 in Fig. 4 (c)). 

Next we take a closer look at the asymptotic behavior at T — > oo. In Fig. 5 we have 
plotted the mean travel distance d defined in the Eq. (Q. For N = 1,2 there are no bulk- 
reptons so the mean travel distances of RD and FM polymers may differ for N > 3 only. 
The calculation reveals that for long RD polymers (N > 5, a 'critical length') the rebound 
effect wins (i.e. d > 0) and the polymer starts traveling backwards while the single particle 
and FM polymers are traveling to the expected negative direction. The rebound effect is 
also present in FM polymers, but it is not strong enough to reverse the drift direction. For 
RD polymers with L > 3 with feasible polymer lengths our Monte Carlo test simulations do 
not display this kind of an anomalous current inversion, suggesting that it may be related 
to spatial discretization and that longer-range interactions (e.g. stiffness) between reptons 
need to be introduced to see such inversions for L > 3. 
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FIG. 4: (Color online) Drift and Peclet numbers of the RD polymers (left column) and FM polymers 
(right column) in the flashing ratchet (L = 3, V mSuX = 1/2) with N = 1 (solid black), N = 3 (dash), 
N = 5 (dot), N = 7 (dash-dot), N = 9 (solid blue). Symmetry parameters x = T q r/T are x = 1/4 
(a), x = 1/2 (b) and x = 3/4 (c). 
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FIG. 5: Mean travel distances of the RD polymers (circles) and FM polymers (squares) in the 
flashing ratchet (L = 3, V max = 1/2) as a function of the polymer length in one time-period at 
asymptotic limit (i.e. the stationary state is reached before the switching). 
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We also note that a similar effect of multi ple current inversions with tightly connected 



Brownian particles (rods) was reported in Ref. 



121 ] . In that work, however, current inversions 



were not found for objects able to vary their length (rotating rods) in the ratchet direction, 
whereas the polymers in our work are able to vary their length between 1...N and still have 
drift inversion. 

The reason for the stronger rebound effect of the RD polymer is caused by the strong ten- 
dency to enter (possibly deformed) U-shaped configurations because of the strict reptation 
rule. After the potential is turned off, this shape unwinds and causes the drift. This also 



happens with time-dependent fields 



221 ] . Since FM polymers lack the reptation rule, there 



is not as much variation in their shape as RD polymers have, thus resulting in a weaker 
rebound effect. 



C. Traveling ratchet potential 

Let now T = T\ + Ti for the mean life-times T\ and T2 of the potentials depicted in the 
right column of Fig. 2 and define the symmetry parameter x = Ti/T. A Similar drift and 



diffusion behavior as previously reported in Ref. |17l . [34| for a single particle is expected. 
In Fig. 6, we show v as a function of x with three different T's (Figs. 6 (a)-(c)): T — > 
0, In (T) = 3 and In (T) = 7. The behavior for the single particle is as expected; the drift 
is antisymmetric with respect to x = 1/2 and goes to zero at x — 0, 1/2, 1. With longer 
polymers the drift changes sign non-trivially for large T's (Fig. 6 (c)) for both polymer types. 
This result is unexpected. An example of the behavior of the diffusion coefficient is shown 
in Fig. 6 (d) for N = 9 and different T's. D c r always reaches its maximum at x = 1/2 and 
decreases as the system goes to a static potential state at x = and 1. The similarity of 
Figs. 3 (a) and 6 (a) is caused by the fact that, as it can be easily seen from Eq. ([H]) for 
T — > 0, the traveling potential creates a similar effective rate structure as the non-symmetric 
flashing ratchet. 

Next we fix x = 1/4 and examine the T dependence in detail. In Fig. 7, we have plotted 
v and the Peclet number for In (T) = —4. ..7. 5. As N > 2, drift inversions can be seen around 
In (T) w 2 for both polymer types. As before, the single particle remains the fastest for small 
T, but eventually the drift curves begin to intersect as T gets larger and the single particle 
is not always the fastest (see e.g. the N = 3 FM polymer in Fig. 7 (a), right column). The 
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FIG. 6: (Color online) Drift and diffusion of the RD- and FM polymers in the traveling potentials 
(L = 3, Vmax = 1/2). (a-c): drift as a function of the symmetry parameter x = T\/T with T — ► 
(a), T = exp(3) (b) and T = exp (7) (c) with N = 1 (solid black), JV = 3 (dash), N = 5 (dot), 
N = 7 (dash-dot), N = 9 (solid blue), (d): the effective diffusion coefficient for N = 9 with T —* 
(solid), T = exp (3) (dash) and T = exp (7) (dot). 

behavior of the Peclet number is as before: Longer polymers have more coherent transport, 
excluding the possible drift inversion points and their neighborhood. With small values of 
T, the Peclet number is also the same for both polymer types, but because of unequal drifts 
for moderate and large values of T (In (T) > 0), also differences exist. 

The insets of Fig. 7 show the drift as a function of iV = 1...11 in detail. We have chosen 
In (T) = 2.85 for RD polymers and In (T) = 2.05 for FM polymers. With these choices, the 
drift inversion occurs between N = 6 and 7 for both models. In the insets of Fig. 7 (b) we 
have plotted the overall drift minimal values in the interval In (T) = —4. ..7. 5 as a function of 
N. The distinction between the polymer types is very clear. FM polymers drift increasingly 
fast backwards whereas RD polymers eventually stop moving as N gets larger. The drift 
inversion of the RD polymers iV > 10 would require a smaller fixed x. 

The magnitude of the drift, typically between 10~ 5 and 10~ 3 , is comparable with the drift 
caused by a flashing ratchet. The Peclet number values of the polymer motion remain small 
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FIG. 7: (Color online) Drift and Peclet numbers of the RD polymers (left column) and FM polymers 
(right column) in the traveling ratchets (L = 3, V mSuX = 1/2) as a function of the mean time-period 
T with the symmetry parameter x = T\/T and N = 1 (solid black), N = 3 (dash), N = 5 (dot), 
N = 7 (dash-dot) and N = 9 (solid blue). For the left inset of (a) In (T) = 2.85 and In (T) = 2.05 
for the right inset. 

(<C 0.1) for both potential types, indicating very low coherence of transport. 



IV. CONCLUSIONS 



We studied the ratchet effect with discrete polymer models in time-dependent potentials 
using the master equation approach. We gave general equations for calculating the effective 
diffusion coefficient and drift in time-dependent periodic systems. Using these equations, 
we performed calculations in the flashing and traveling ratchet potentials for short discrete 
polymers with the Rubinstein-Duke model and a relaxed version of this model allowing tube 
breaking. We found complex dynamics that results from the non-pointlike structure of the 
polymers by the coupling between the potential and polymer internal states. By varying 
the potential switching rates, we found non-trivial inversions of the polymer drift direction, 
which cannot occur with simple pointlike non-interacting particles. We also found that the 
Peclet number grows as the polymer gets longer and is largely independent of the polymer 
type thus allowing more coherent transport for longer polymers. The overall polymer dy- 
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namics in ratchet potentials was found to be very model specific. The discretization of the 
problem in this work may be far from many real- world applications but, nevertheless, since 
our model catches the essential characteristics of the Brownian motor system, we expect 
that similar properties could be found in the nano-scale objects that can be described with 
discrete states instead, such as molecular motors with internal structure. Drift inversions 
are especially interesting since they facilitate more efficient separation methods of molecules. 
The next step would be to consider larger L and N and the differences between deterministic 



and stochastic switching [35]. 



APPENDIX A: OPERATORS IN H 

The explicit definitions of the operators in Eq. (1) are 

Md) = {R(l) + L(l)}h ,i,i - R(l)~a\ tl - L(/)6 f M 

+ L(l)n Axl - L(l)a 1:l + R{l)n B ,i,i ~ R(l)h,i 

B y>l = {R(l + f(N - 1, y)) + L(l + f(N - 1, y))}n 0;N ^ y>l 

- R(l + f(N - 1, y))oLi, v> , " L(l + f(N - 1, 2/))&Li^ 
+ L(l + f(N - 1, y))nA,N-i, v ,i ~ Hi + f(N - 1, y)) 

a N-l,y,l 

+ R(l + f(N - 1, j/))n fllJV _i, V) , - R(l + f(N - 1, y))b N . ltVtl 
M itVt i = R(l + f(i,y))(n A ,i, y ,in 0i i + i tyj i + no t i >y! in B ,i+i, y ,i - ai, y ,i(4+i, y ,i ~ b i, y ,fii+hv,i) 

+ L(l + f(i,y))(n ,i,y,l n A,i+l,y,l + n B,i,y,l n O,i+l,y,l a i,y,l a i+l,y,l bi : yjb i+lyl ) 

+ QR(l + f(i, y)){n A ,i, y ,inB,i+i, y ,i + no.i.y.jno.i+i,^ - a^bi+x^i - 4,y,j a i+i,j/,z) 

,i,y,l n A,i+l,y,l + ^0,j,j/,« n 0,i+l,y,/ bi^jdi+i^j a \ y y ^y , 
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where f2 = for RD polymers and 1 for FM polymers, and 

a hl = c+ci! a{ tl = cja\ 

k,i = cf&i b\ :l = c+b{ 

n z ,i,i = nin z ,i 

/i-i \ 

x i,y,l = n l I Y\ n 9{y,j)J I x * 
n z,i,y,l = ni I J^J I ^z,i 

with a: G {a, 6, a*, 6^}, z G {A, 0, 5}. The function i) G {A, 0, 5} gives the state of the 
ith bond in the configuration y, and the function / 

i 

f(i,y) = 1 $2(S!y\n Ati -nB,i\*y), l<t<N-l 

3=1 

gives the position of the repton i + 1 in marker-centered coordinates. The detailed forms of 
the functions g and / depend on the selection of the state basis. 
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